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Abstract 

In a recent paper [A. Moujahid, A. d'Anjou, F. J. Torrealdea and F. Torrealdea, Phys. Rev. E 
. 83, 031912 (2011)], the authors have calculated the energy consumed in firing neurons by using 

the Hodgkin-Huxley (HH) model. The energy consumption rate adopted for the HH model yields 

Ch ' 

^ ' a negative energy consumption meaning an energy transfer from an HH neuron to a source which 

is physically strange, although they have interpreted it as a biochemical energy cost. I propose an 
', alternative expression for the power consumption which leads to a positive energy consumed in an 

HH neuron, presenting some model calculations which are compared to those in their paper. 

O ' PACS numbers: 87.19.11, 87.19.1g,87.19.1y,87.18.Sn 
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It has been controversial how to experimentally and theoretically evaluate the ener; 
consumption in firing neurons. Recently Moujahid, d'Anjou, Torrealdea and Torrealdea [l| 
have proposed such a theoretical method by using the Ho dgkin- Huxley (HH) model given 

by 

CV = -/Na - /k - /l + /, (1) 

m = -{am + bm)m + a^, (2) 
h = -{ah + bh)h + ah, (3) 
n = -{an + bn)n + an- (4) 

Here V denotes the membrane potential in mV, / stands for the total current in /xA/cm^, 
the membrane capacitance is C = 1 /iF/cm^, and m, h and n are dimensionless gating 
variables of Na, K and leakage (L) channels, respectively. Currents in respective channels 
are expressed by 

In^ = 9Ka.m^h{V -V^^), (5) 
Ik = 9Kn\V-VK), (6) 
lL = gL{V-Vi^), (7) 

where reversal potentials of Na, K and L channels are l^Na = 50 mV, Vk = —77 mV and 
Vl = —54.5 mV, respectively, in the conventional absolute unit [2I, Q], and the maximum 
values of corresponding conductances are (/Na = 120 mS/cm^, (/k = 36 mS/cm^ and ol = 0.3 
mS/cm^ 0,0]. Coefficients of a.m and bm et. al. are given in the absolute unit by 

am = 0.1(V^ + 40)/[l-e-(^+^°)/i°], 

_ 4„-(F+65)/18 

ah = 0.07e-(^+65)/2o^ 

bh = l/[l + e-(^+3^)/^°], 

an = 0.01(r + 55)/[l-e-(^+^^)/i°], 

6„ = 0.125 e"(^+6^)/«°. 



In Ref.pJ the total energy H is expressed by 
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H = -CE' + i/Na + H^ + Hj^, (8) 



where the first term expresses the energy stored in the capacitor and Hi denotes the energy 
in the channel i (= Na, K and L). Based on a biochemical consideration, the authors in 
Ref. 1] have derived a derivative of H with respect to time given by |l| 

H = P' = CEE + /Na^Na + /r^K + /l^L, (9) 

where E {= V — l^es) and Ei {= Vi — l^es) are action potential and reversal potentials, 
respectively, in the reduced unit with V^es (= —65 mV) [3]. The four terms in Eq. ([9]) stand 
for energy consumption rates in respective channels. By using Eq. ([H]), the authors in Ref. |l| 
have studied the energy consumption rate against firing rate in an HH neuron. Equation ([9]), 
however, yields a negative energy consumption [Fig. 2(b) in Ref. for which alternative 
methods are studied as will be explained in the following. 

It is well known that the HH equation given by Eq. ([T]) expresses an electric circuit (see 



Fig. 1 in Ref. [2]) consisting of capacitor C and three Na, K and L ion channels, which are 
connected in parallel. The channel i (= Na, K and L) includes a resistor Ri and a battery 
Vi, through which a current Jj flows. The total current / flows when these components are 
connected to a source battery V . Applying Kirchhoff's law to the circuit, we obtain 

/ = Cy + ^/i, (10) 

i 

liRi = V -Vi (i = Na, K and L). (11) 

Equation ( ITOll is nothing but Eq. ([T]). Total consumed energy rate (power) in the circuit is 
given by 

p = cVV + J2iPJ^ + PR^, (12) 

i 

= CVV + J2Pi, (13) 

i 

with 

Pji = lfR, = h{V-Vi}, (14) 
Pm = Wi, (15) 

where IfRi and liVi signify contributions from Joule heat and reversal potential, respectively, 
in the channel i. 
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Now we consider three methods A, B and C for a calculation of the power consumed in 
an HH neuron, depending on which contributions are taken into account, 



= cVV + J^PRi (method A), (16) 

i 

PB = CVV + J2Pji (method B), (17) 

i 

Pc = CVV + 5^(Pffi + Pji) = VI (method C). (18) 

i 

Methods A, B and C take into account contributions from the reversal potential {Pm), Joule 
heat (Pji), and the reversal potential plus Joule heat (Prj + Pji), respectively, besides that 
from the capacitor (CVV). The VI term in Eq. ( fTSi) expresses a total power supplied from 
a source. By a simple calculation, the total power adopted in Ref. ij] becomes [Eq. (^] 

P' = Pa- Kes /, (19) 

which does not include a contribution from Joule heat. 

The inset of Fig. 1(a) shows the time course of an action potential for an external current 
of / = 6.9 yuA/cm^, for which a neuron fires with a period of 17.36 ms [l]. Properties of the 
consumed power depend on which method among the three is adopted. Figure 1(a), 1(b) 
and 1(c) show time courses of power consumption in four channels and total power for a 
single firing calculated in the methods A, B and C, respectively. We note that total power 
Pa in Fig. 1(a) is negative, Pb in Fig. 1(b) is positive, and Pq in Fig. 1(c) is oscillating 
because Pc = IV. 

From the total power P, we may evaluate the mean total power P, which is given by 
Pa < 0, Pb > and Pc < for methods A, B and C, respectively, a bar denoting an average 
over a period. Absolute mean powers \Pk\ (k = A, B and C) calculated in the methods A, B 
and C are plotted by solid, dashed and chain curves, respectively, as a function of external 
current / in Fig. 2, whose inset shows the / dependence of firing frequency / of an HH 
neuron. The mean power of firing state in the method A is about 10000 ~ 15000 nJ/s for 
/ = 7 ~ 30 yuA/cm^ while that of the quiescent state without firings is much smaller (300 
~ 900 nJ/s). I Pa I shown by the solid curve is similar to the result presented in Fig. 3(b) of 
Ref.|l|. The mean power in the method B is almost the same as that in the method A, 
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FIG. 1: (Color online) Time courses of powers consumed in a capacitor (C), Na, K and L channels 
and total powers calculated in (a) the method A, (b) method B and (c) method C with 1 = 6.9 
/iA/cm^, a total power in (c) being multiplied by a factor of ten. The inset in (a) shows the time 
course of an action potential V . 



which arises from the relation, 



(21) 



because VI ^ | Xli-^ral (Fig- !)• The mean power in the method C is smaller than those 
in the methods A and B, and it is almost linearly increased with increasing J. The positive 
mean power implies that energy is supplied from an external source to a neuron while the 
negative power means the opposite in which a neuron plays a role of a power generator. The 
authors in Ref. jll have interpreted an obtained negative energy consumption rate [P' < 0) 



i6ooa 



12000- 
% 8000- 

1^ 

4000h 



B 



C(xio) ./■ga 



^ 60- 

N 

^ 40- 
20h 



^ 10 20 2 30 
I I {\xhlcva ) 



10 



20 



30 



I (|iA/cm^) 



FIG. 2: (Color online) Absolute mean total power |P| as a function of input current I calculated 
in the method A (solid curve), method B (dashed curve) and method C (chain curve), result of 
the method C being multiplied by a factor of ten. The / dependence of firing rate / is shown in 
the inset where a firing occurs at / > 6.2 /iA/cm^. 

as a loss in electrochemical energy which is expected to correspond to consumed power 
because it has to be recharged by adenosine triphosphatase (ATPase) activity. Equation 
( pTl) shows that the energy consumption assumed in the method A (also in Ref. [l|) may be 
related to Joule heat consumed in conductors of ion channels in the model B. 

The authors in Ref. [l| have extended their analysis to two-electrically coupled HH neu- 
rons and also to a collection of uncoupled HH neurons, for which expressions for energy 
consumption rates similar to Eq. ([9]) have been derived {e.g. Eq. (8) in Ref. {l| for two-HH 
neurons). They have the same deficit as in the case of a single HH neuron having been 
mentioned above. 

To summarize, the calculation in Ref. [l| based on Eq. ^ yields a negative mean power 
consumption in an HH neuron, which is physically inappropriate, although its biochemical 
interpretation is reasonable. We have proposed an alternative method, considering three 
methods A, B and C given by Eqs. f|T6l) - f|T8|) . The method B which takes into account a 
contribution from Joule heat consumed in conductors, leads to a positive mean power in the 
HH model. Energies calculated in the reduced unit adopted in Ref. 1| are different from 

n p _ M 

those in the absolute unit [3!]. For example, the mean value of EI in the reduced unit is 
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positive Q while its counterpart VI in the absolute unit is negative. By using the method 
B, we may study nonequilibrium properties of an HH neuron from a viewpoint adopted in 
Ref. where an analogy between an electric RC circuit and a damped Brownian particle is 
employed. 
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